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Abstract 

LOFAR is a new form of radio telescope which can detect radio emission 
from air showers induced by very high-energy cosmic rays. It can also look 
for radio emission from particle cascades on the Moon induced by ultra high- 
energy cosmic rays or neutrinos. To complement the radio detection, we are 
setting up a small particle detector array LORA (LOfar Radboud Air shower 
array) within an area of ~ 300 m diameter in the LOFAR core. It will 
help in triggering and confirming the radio detection of air showers with the 
LOFAR antennas. In this paper, we present a short overview about LORA 
and discuss its current status. 

1 Introduction 

LOFAR (the LOw Frequency ARray) is a new kind of radio telescope for astro- 
nomical observations in the frequency range of ~ (10 — 240) MHz with high 
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Figure 1 : Detectable CR energy ranges for LOFAR. LOFAR VHECR refers to 
the detection using air showers and UHEP to the detection using the Moon. Also 
shown are the energy ranges for the LOPES and the Auger radio experiments. 

sensitivity and high spatial resolution (http://www.lofar.orgy It uses a large num- 
ber of simple dipole antennas instead of the traditional big parabolic dishes. It 
consists of 40 stations in the Netherlands, 5 in Germany and one each in Great 
Britain, France and Sweden covering a total area of more than 1000 km in diame- 
ter. 

Though primarily design as an astronomical telescope, LOFAR can also be 
used for the detection of very high-energy cosmic rays (CRs) in the interesting 
energy region above 10^^ eV where the transition of galactic to extra-galactic CRs 
is expected (Horandel et al. 2009, Homeffer et al. 2010). This will be done 
by looking at extensive air showers which are essentially cascades of energetic 
secondary particles produce by the interaction of CR primaries with the nuclei 
present in the atmosphere. A large fraction of these secondaries are electrons and 
positrons which produce radio synchrotron emission in the presence of the Earth's 
magnetic field (Falcke et al. 2005). Due to coherence effects, this emission can 
give strong signals on the ground in the frequency range of ~ (10 — 80) MHz 
detectable by the LOFAR low band antennas. 

In addition, LOFAR can also look for the highest energy CRs or neutrinos 
above around 10^^ eV (Scholten 2010). This will be carried out by detecting co- 
herent radio Cherenkov emission from particle cascades at the Moon induced by 
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those CRs or neutrinos. At frequencies of ~ (100 — 200) MHz which is within the 
range of the LOFAR high band antennas, the angular spread of this emission be- 
comes wider leading to more detectable signals at the Earth. The different energy 
ranges detectable by LOFAR are shown in Fig. 1 where the VHECR (Very High 
Energy CRs) refers to the detection using air showers and the UHEP (Ultra High 
Energy Particles) refers to that using the Moon. Fig. 1 also shows the energy 
ranges for the LOPES (LOFAR Prototype Station) experiment which is located 
at the KASCADE-Grande experimental site (Apel et al. 2010) and the currently 
building AERA (Auger Engineering Radio Array) experiment at the Pierre Auger 
Observatory (van den berg et al. 2009). 

One important goal of the LOFAR CR experiment is to push the radio detec- 
tion technique towards an independent way of detecting very high energy CRs. 
By detecting the radio signals with better sensitivity and better spatial resolution 
in a wider frequency range, we strongly believe that LOFAR will provide better 
understanding of the measured signals, their emission mechanisms and their re- 
lations with the air shower parameters, thereby leading to better estimates of the 
properties of the primary particle. Compared to the LOPES (30 antennas) and the 
CODALEMA (24 antennas) experiments (Ardrouin et al. 2005), LOFAR has 18 
stations in its core (an area of 2 x 3 km^) with each station consisting of 96 low 
band antennas and 48 high band antennas. However, at this stage, it is still quite 
early for radio detection experiments to do a stand-alone study on CRs. There- 
fore, we are also setting up a small particle detector array called LORA (LOfar 
Radboud Air shower array). Its main role will be to trigger and confirm the radio 
detection of air showers with the LOFAR antennas. It will also help in the recon- 
struction of several important air shower parameters like the primary energy, the 
shower core location, the arrival direction, the lateral density distribution etc. 

2 LORA set-up 

LORA is setting up inside the LOFAR core within an area of ~ 300 m diameter. 
It consists of 5 stations with 4 particle detectors each, placed at a separation of ~ 
(50 — 100) m between them. It is expected to detect CRs with energies greater than 
~ 10^^ eV at an event rate of around once every few minutes. Detailed simulations 
about LORA as well as combined studies including the LOFAR radio data will be 
carried out as soon as possible. The schematic layout of the LORA detectors 
arrangement is shown in Fig. 2 where the red stars represent the positions of 
the detectors and the circles denote the positions of the LOFAR antennas fields. 
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Figure 2: Layout of the LOfar Radboud Air shower array (LORA) in the LOFAR 
core. The red stars denote the positions of the particle detectors and the circles 
denote the positions of the antenna fields. The detectors are contained within an 
area of ~ 300 m diameter with relative spacings of ~ (50 — 100) m. 

Data collection in individual stations are controlled locally by station computers 
which are then controlled by a master computer where the overall data processing 
is done. 

2.1 Detectors 

The detectors for LORA are taken over from the KASCADE experiment (Antoni 
et al. 2003). Each detector is of (98 x 125) cm in size and consists of two slabs of 
3 cm thick plastic scintillators. When a charged particle passes through a scintil- 
lator slab, the light which are produced are collected into a photomultiplier tube 
(PMT) through a wavelength shifter bar. The electrical signals from the PMTs 
(one PMT for each slab) are then combined and fed into the electronics where 
they are converted into digital signals. 
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Figure 3: Example of a signal trace when a charged particle passes through one 
of our detectors. 

2.2 Electronics 

The electronics was developed originally for the HISPARC experiment (http://www.hisparc.nr). 
An electronics unit (hereafter referred to as HISPARC unit) can handle two chan- 
nels (one channel corresponds to one detector) and there are two units per station 
(a master and a slave). The master sets the event trigger condition and also pro- 
vides the time stamp with a GPS receiver along with a 200 MHz clock counter. 
Data in each channel are handled by two 12-bit ADCs (so there are 4 ADCs per 
HISPARC unit) which can measure voltages in the range of (0 — 2) V. Event data 
(signal) in a channel are sampled with a time resolution of 2.5 ns and are stored in 
a total time window of 10 /iS. Fig. 3 shows a typical signal of a charged particle 
passing through one of our detectors. Data from the electronics are sent to the 
station computer through USB. 

Each HISPARC unit has a FPGA circuit built into it. This gives an observer 
control over several parameters required by the electronic s as well as the detectors. 
For instance, one can set the trigger thresholds, high voltages, trigger condition, 
coincidence time window etc. at any time during the observation through a user 
software running at the station computer. 
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2.3 Data acquisition 

AH the PMTs are gain calibrated by adjusting their voltages to match the response 
for single muons. When the signals in the two HISPARC units in a station sat- 
isfy some minimum trigger condition, the signals are sent to the station computer 
which then sends them to the master computer. The master computer combines 
the signals from all stations and checks for good air shower events. When an air 
shower candidate is found, a first-level online analysis is performed to calculate 
shower parameters like the arrival direction, shower core position, event size etc., 
and a trigger information will be sent to the LOFAR antenna to dump the cor- 
responding radio data for the air shower. The master computer takes less than 
10 ms to process an event. This is necessary because the antennas have to dump 
their radio data corresponding to a particular event before they are overwritten in 
a memory ring-buffer, known as Transient Buffer Board (TBB). 

The DAQ software for LORA is developed in C/C++ language for Linux based 
operating systems. The software (particularly the online monitoring tool) uses 
several features of the ROOT package (http://root.cem.ch). On the monitoring 
panel, important information is displayed which is useful to monitor the perfor- 
mance of the system electronics and of the detectors during observations. 

The final data are stored in ROOT format and they consist of four kinds of 
data. The first are the event data which are generated whenever an air shower 
event is detected. The data contain the event time stamp for each station and the 
signal trace in ADC counts for each detector. The second kind of data stores the 
so-called one second messages from the HISPARC unit. This data is generated 
every second by the master device and it contains information about the number 
of times the analog signal went over the threshold in the last second for each of 
the four channels. It also has important timing information which can be used for 
calculating an event time stamp with nanosecond accuracy. The third data kind 
comprises several control parameters for the observation run. This data is stored 
every interval of time fixed by the observer at the start of the run. The fourth kind 
of data contains information about the noise level in each channel averaged over 
some fixed interval of time. 

The DAQ software and a preliminary data analysis software have been tested 
successfully on a LORA prototype installed at Radboud University (RU) Ni- 
jmegen. The results from the prototype are presented in the following section. 
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Figure 4: Reconstructed air shower core distribution for 5463 events collected 
with the LORA prototype. 

3 Results from the LORA prototype 

The LORA prototype that we set up to test the electronics and the DAQ software 
consists of 4 detectors in a (20 x 80) m arrangement. With that, we collected 
more than 5000 air shower events in a total observation time of ~ 236 hrs (~ 10 
days) during April 2010. This corresponds to an event rate of around once every 
3 minutes. We have performed a preliminary analysis of the data to reconstruct 
the shower core position on the ground, the arrival direction and the energy of 
the primary. The primary energy is calculated from the total number of charge 
particles on the ground using the parameterization given in Horandel 2007. This 
simple estimate gives the energy threshold of the prototype to be less than 10^^ eV. 
The shower core position is calculated using the center of gravity method and the 
primary arrival direction is calculated using the relative pulse timing informations 
in the detectors assuming that the shower front arrives in a plane on the ground. 
The normal to the shower plane gives the arrival direction and it is characterized 
by two angles: the zenith angle 9 measured from the vertical direction and the 
azimuthal angle (f) measured in the horizontal plane. Though detailed simulation 
is still yet to be performed, the energy resolution of LORA is expected to be 
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Figure 5: Arrival direction distribution {9, (p) for the 5463 events collected with 
the LORA prototype. 

< 30% and the angular resolution to be < 1°. 

In Figs. 4 and 5, we show the distributions of the reconstructed core positions 
and the arrival directions (9, 0) respectively. In Fig. 6, we plot the 9 (top) and the 
(p (bottom) distributions separately. The 9 distribution is fitted using the following 
function, 

^ = a sin6' cos''^. (1) 
d9 

The fit parameters are found to be a = 636 ± 14.9 and b = 6.35 ± 0.14. For the cp 
distribution, we fit it to a constant as 

dN 

Eq. (2) is expected if the arrival direction of the CRs are isotropic. However, in 
our case we see a sinusoidal distribution with two dips one at ~ 90° and the other 
at ~ 270° which are due to the narrow arrangement of our detector set-up. 

It should be mentioned that at this stage, we do not aim to perform a detailed 
analysis and derive a final set of fit parameters which determine important proper- 
ties of the CRs. As already mentioned, our primary aim was to test our electronics 
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and the DAQ software. The preUminary results from the prototype are in general 
consistent with what we expect from CR observations with air shower arrays. 

4 Current status and future plans 

The electronics and the DAQ software for LORA have been tested successfully. 
The installation of all the LORA stations in the field is already completed and 
the set-up is currently under testing. On the other hand, the LOFAR CR pipeline 
software is under development and we are expecting to start the simultaneous ob- 
servation of CRs with the radio antennas before summer 20 11. 
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Figure 6: Zenith 9 (top) and azimuthal angles (bottom) distribution for the mea- 
sured showers. The 6 distribution is fitted with a function given by Eq. (1). The 
distribution is fitted with a constant. 
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